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I .  INTRODUCTION 


As  part  of  its  research  effort  on  impact  acceleration  injury 
prevention,  the  Naval  Biodynamics  Laboratory  (NBDL)  has  collected 
a  fairly  extensive  data  base  of  93  -G^  accelerator  runs  using  rhesus 
monkeys  as  subjects.  Desmatics  has  used  these  data  to  formulate  sev¬ 
eral  statistical  models  which  predict  the  probability  of  head/neck 
impact  acceleration  injury.  The  primary  goal  of  this  line  of  research 
is  the  establishment  of  human  tolerance  limits  to  -G^  acceleration. 
These  tolerance  limits  will,  in  turn,  provide  essential  input  to  the 
process  of  evaluating  proposed  protection  systems  for  naval  aircraft. 

Any  statistical  injury  prediction  model  developed  from  data  for 
a  single  primate  species  will  be  inadequate  for  the  purpose  of  pre¬ 
dicting  injury  for  other  primate  species.  However,  if  data  from  two 
or  more  species  are  used  to  develop  a  composite  model  which  accurately 
fits  the  experimentally-obtained  injury  data  from  each  species,  empir¬ 
ical  information  about  species-specific  effects  is  obtained.  This 
information  should  then  be  useful  in  efforts  to  extrapolate  the  injury 
prediction  model  to  humans. 

The  only  previous  primate  experiments  using  -Gx  indirect  impact 
acceleration  were  carried  out  by  Clarke  et.  al.  [3-5]  with  baboons. 

They  investigated  tolerance  to  abrupt  linear  deceleration  using  various 
restraint  systems.  Only  their  last  series  of  experiments,  as  described 
in  [3],  is  of  interest  here.  Those  experiments  used  the  Air  Force 
shoulder  harness-lap  belt  restraint  and  the  primary  cause  of  death  was 
lower  brain  stem  or  cervical  spinal  cord  trauma.  Since  these  injuries 


are  conyarable  to  those  found  by  NBDL  scientists  for  rhesus  monkeys, 
it  seems  reasonable  to  assume  that  the  injury  mechanism  is  the  same. 
In  this  report  we  develop  a  statistical  injury  prediction  model  using 
the  data  for  both  rhesus  monkeys  and  baboons.  This  model  illustrates 
both  the  similarities  and  the  differences  in  the  responses  of  the  two 
species  to  -Gjj  acceleration. 


II.  METHODS 


Several  previous  Desmatics  technical  reports  have  discussed  the 
use  of  a  logistic  function  in  the  development  of  impact  acceleration 
injury  prediction  models.  These  functions  are  of  the  form: 


k 

P(x)  -  {1  +  exp[-(8n  +  E  B.x  )]} 

i«l  1  1 


-1 


where : 

x  m  (Xj,  x^,  ....  x^)  denotes  the  set  of  predictor  variables, 

(8q,  ....  8^)  denotes  a  set  of  parameter  values, 

and  P(x)  denotes  the  true  probability  of  injury  corresponding  to  x. 

Desmatics  has  considered  several  variables  as  possible  predictors 
of  injury.  Empirical  data  have  been  used  to  fit  the  models  and  evaluate 
their  ability  to  accurately  predict  injury.  The  most  recent  Desmatics 
technical  report  [1]  in  this  area  considers  a  model  based  on  peak  sled 
acceleration  and  the  initial  yaw  angle  of  the  head.  That  model  is 
compared  to  earlier  models  and  shown  to  be  a  better  predictor  of  injury. 
Therefore,  in  this  investigation  we  use  peak  sled  acceleration  and  ini¬ 
tial  yaw  to  predict  injury  for  rhesus  monkeys.  These  variables  have 
been  measured  for  21  high-level  accelerator  runs  and  are  given  in 
Table  1. 

Clarke  et.al.  conducted  their  tests  on  the  Daisy  Decelerator  [2] 
with  stopping  distances  varying  from  0.5  to  3.5  ft  at  6  in  increments. 
The  peak  sled  accelerations  for  those  tests  are  given  in  Table  2.  Prior 


Run 

Number 

Subj ect 
Number 

In j  ury 
Code 

LX0660 

A03146 

0 

LX0661 

A03146 

1 

LX1359 

A04099 

0 

LX 1360 

A04099 

1 

LX1362 

A03935 

0 

LX1363 

A03935 

1 

LX1365 

A03921 

1 

LX1893 

A03924 

0 

LX 1894 

A03933 

0 

LX1895 

A03951 

1 

LX1896 

A03946 

1 

LX1905 

A04101 

1 

LX3015 

AR0764 

0 

LX3016 

AR0764 

1 

LX3026 

A03923 

1 

LX3031 

AR4115 

1 

LX3032 

AR3936 

0 

LX3033 

AR3936 

1 

LX3188 

AR8863 

0 

LX3192 

AR8866 

1 

LX3709 

AR8790 

1 

Peak  Sled 
Acceleration  (g) 

Initial 
Yaw  (deg) 

107 

0 

158 

0 

110 

-30 

128 

90 

108 

30 

123 

0 

109 

90 

110 

30 

108 

-30 

131 

-90 

131 

0 

126 

0 

97 

90 

124 

90 

127 

-60 

127 

0 

125 

0 

163 

0 

104 

-30 

105 

60 

88 

60 

Table  1.  High-Level  Rhesus  Accelerator  Runs. 


Injury  Code»l  Denotes  Fatal  Run 


Injury  Peak  Sled  Injury  Peak  Sled  Injury 

Code  Deceleration  (g)  Code  Deceleration  (a)  Code 


0.5 

ft 

Stopping  Distance 

0 

21.15 

0 

44.14 

0 

29.02 

0 

61.16 

0 

35.03 

0 

67.27 

0 

42.14 

0 

81.99 

1.0 

ft 

Stopping  Distance 

0 

11.34 

0 

25.38 

0 

16.96 

0 

34.17 

0 

24.86 

1.5 

ft 

Stopping  Distance 

0 

12.96 

0 

32.13 

0 

24.80 

0 

40.32 

2.0 

ft 

Stopping  Distance 

14  tests 

6.50 

0 

69.36 

no  injury 

to 

26.19* 

0 

82.46 

0 

35.89 

0 

87.01 

0 

43.87 

0 

87.79 

O 

44.47 

1 

92.86 

0 

48.33 

1 

97.14 

0 

56.63 

1 

97.42 

0 

60.80 

0 

102.1 

0 

67.52 

1 

105.1 

2.5 

ft 

Stopping  Distance 

0 

11.43 

0 

32.19 

0 

21.34 

0 

37.43 

3.0 

ft 

Stopping  Distance 

0 

8.08 

0 

26.04 

0 

20.00 

0 

38.82 

3.5 

ft  Stopping  Distance 

0 

9.70 

0 

60.63 

1 

0 

18.97 

0 

74.61 

1 

0 

26.11 

1 

89.76** 

1 

0 

37.79 

0 

90.36 

1 

0 

_  50.52 

0 

97.42 

1  *Not  used  In  logistic  estimation  sines  exact  accelerations  unknown. 

[  **lep  belt  severed  during  lapset.  Not  used  in  logistic  sstinatlon. 

i  Table  2:  Baboon  Decelerator  Runs 

i 

i 


Peak  Sled 
Deceleration  (g) 


98.88 

113.6 
118.4 

133.6 


45.74 

54.71 


52.20 

61.20 


107.3 
109.9 

113.4 

113.6 
118.0 
118.2 
121.8 

130.6 


53.01 

59.44 


52.07 

62.22 


99.12 

100.80 

111.36 

112.10 


Injury  Code*l  Denotes  Fatal  Run 


to  experimentation,  each  animal's  head  was  fixed  in  place  with 
masking  tape,  eliminating  any  dependence  of  injury  on  initial  head 
position.  We  therefore  consider  only  stopping  distance  and  peak 
sled  acceleration  as  possible  predictions  of  injury  for  the  baboon 
experiments . 


III.  RESULTS  AND  DISCUSSION 


Clarke  et.al.  fit  probit  models  to  their  data  and  conclude  that 
impact  acceleration  injury  depends  solely  on  peak  sled  acceleration 
for  the  range  of  impact  parameters  used  in  their  study.  They  base 
this  conclusion  on  the  fact  that  the  calculated  LD5QS  (the  sled  accel¬ 
erations  where  injury  is  predicted  for  50%  of  the  subjects)  vary  little 
over  the  range  of  stopping  distances.  Here  we  conduct  a  formal  signi¬ 
ficance  test  of  the  hypothesis  that  the  probability  of  injury  does 
not  depend  on  stopping  distance. 

If  stopping  distance  has  no  effect  on  the  probability  of  injury 
then  a  single  logistic  model  based  on  all  of  the  data  should  fit  as 
well  as  a  model  which  employs  a  separate  logistic  function  for  each 
stopping  distance.  We  therefore  calculate  the  following  quantities: 

Lj  -  -2  log  likelihood  for  the  model  using  only  peak  sled 
acceleration 

L 2  ■  -2  log  likelihood  for  the  model  which  considers  both 
peak  sled  acceleration  and  stopping  distance. 

Under  the  null  hypothesis  that  stopping  distance  has  no  effect, 

L^-L2  has  an  approximate  Chi-square  distribution  with  12  degrees  of 

freedom.  For  this  data  set,  the  calculated  quantities  are: 

Lr^. 3,  L2-13.6,  LrL2-5.7. 

The  90C^  percentile  of  a  Chi-square  distribution  with  12  degrees  of 
freedom  is  18.55.  Since  the  calculated  statistic  is  less  than  this 


value,  there  is  no  evidence  that  stopping  distance  affects  the  prob¬ 
ability  of  injury.  We  therefore  use  only  peak  sled  acceleration  as  a 


predictor  of  fatality  for  baboons. 

In  order  to  compare  the  probability  of  injury  for  rhesus  monkeys 
and  baboons,  we  fit  the  following  model: 

P(I,A,Y)  -  (I  +  exp[-(e0+61-I+62-A+B3-I-A+64-Y)]}_1 

where : 

I  *  0  if  rhesus  monkey,  I  if  baboon, 

A  -  peak  sled  acceleration, 

Y  ■  initial  yaw  angle  of  the  head, 
and  (8q,  ,  ....  8^)  are  parameters  to  be  estimated. 

This  model  fits  separate  logistic  functions  of  peak  sled  acceleration 
for  the  two  species  and  also  includes  the  effect  of  the  initial  yaw 
angle  of  the  head.  One  question  of  interest  concerning  the  model  is 
whether  the  estimate  of  8^  is  significantly  different  from  zero.  The 
following  quantities  are  calculated: 

Lj  ■  -2  log  likelihood  for  the  model  with  8^  =  0 
and  L2  ■  -2  log  likelihood  for  the  model  with  8^  unrestricted. 

Under  the  null  hypothesis  that  8^  ■  0,  Lj-L2  has  an  approximate  Chi- 
square  distribution  with  1  degree  of  freedom.  The  calculated  values 
are: 

1^-36.85,  L2-36.78,  Lj-I^-O.O?. 

Clearly,  8^  makes  no  important  contribution  to  the  model.  Therefore, 
the  reduced  model  (with  8^*0)  may  be  used  to  describe  the  data.  The 
fitted  model  is  given  by: 

P(I,A,Y)  -  {1  +  exp [-(-22. 204+4. 4293I+. 18202A+.054381Y) ] }_1 
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One  consequence  of  removing  6^  from  the  model  is  that  the  logis¬ 
tic  curves  for  rhesus  monkeys  and  baboons  never  cross.  Furthermore, 
the  differential  effect  of  peak  sled  acceleration  is  the  same  for  the 
two  species.  That  is,  the  increase  in  peak  sled  acceleration  needed 
to  raise  the  probability  of  injury  from  x  to  x+A  is  the  same  for  each 
species,  for  any  value  of  x.  This  agrees  with  the  conclusions  reached 
by  Smith  [6]  on  theoretical  grounds. 

The  fitted  logistic  response  curves  are  plotted  in  Figure  1  for 
each  species  and  various  values  of  initial  yaw.  For  ease  of  reference, 
the  LD^qS  have  been  marked  for  each  curve.  It  can  be  seen  from  the 
plot  that  for  an  initial  yaw  angle  of  zero  degrees,  the  LD,-q  for 
rhesus  monkeys  is  approximately  24  g  higher  than  that  for  baboons. 

In  addition,  the  LD,.q  for  rhesus  monkeys  is  seen  to  decrease  about 
9  g  for  every  30  degrees  of  initial  yaw.  Recall  that  the  effect  of 
initial  yaw  on  the  probability  of  injury  for  baboons  could  not  be  de¬ 
termined  because  each  animal's  head  was  taped  prior  to  experimentation. 

It  should  be  noted  that  the  estimated  logistic  response  curve 
given  here  for  rhesus  monkeys  differs  slightly  from  that  reported  in 
[1].  This  difference  results  from  the  use  of  both  rhesus  and  baboon 
data  to  estimate  the  parameters  in  the  current  model.  However,  since 
the  differentail  effect  of  peak  sled  acceleration  is  not  significantly 
different  for  the  two  species,  there  is  very  little  difference  between 
the  two  models.  The  estimated  response  functions  are: 

P^A.Y)  -  {1  +  exp[-(-18.5+.152-A+.0442-Y)]}~1 
and  P2(A,Y)  -  {1  +  exp[-(-22 . 2+. 182 -A+.0544-Y) ] }_1 . 


-9- 


It  is  very  difficult  to  evaluate  the  effect  of  small  changes  in  the 
estimated  parameters  of  a  logistic  model ,  especially  when  several 
parameters  are  changed  simultaneously.  Figure  2  shows  plots  of  both 
estimated  response  curves  as  a  function  of  peak  sled  acceleration  with 
the  initial  yaw  angle  set  to  zero.  It  can  be  seen  from  the  plot  that 
there  is  little  difference  between  the  two  response  functions.  The 
maximum  difference  between  them  is  less  than  .05.  Furthermore,  the 
estimated  effect  of  the  initial  yaw  angle  of  the  head  is  similar  for 

/N 

the  two  models.  For  P^(A,Y),  the  estimated  LD,.q  decreases  by  8.72  g 
for  every  30°  of  initial  yaw  while  for  P^CA.Y),  the  decrease  is  8.98  g. 
Clearly,  the  differences  between  the  two  models  are  of  no  practical 
significance. 


IV.  SUMMARY 


In  this  report  we  have  compared  data  from  NBDL  -Gx  experiments 
on  rhesus  monkeys  with  the  results  of  Clarke  et.al.  on  baboons  using 
the  Air  Force  shoulder  harness-lap  belt  restraint  system.  Although 
the  restraint  systems  differ  for  the  two  sets  of  experiments,  the 
types  of  injury,  and  by  implication  the  injury  mechanisms,  are  similar. 

In  an  earlier  Desmatics  technical  report,  peak  sled  acceleration 
and  the  initial  yaw  angle  of  the  head  were  shown  to  be  the  best  pre¬ 
dictors  of  injury  for  rhesus  monkeys.  Those  variables  were  therefore 
used  in  the  model  developed  here.  For  the  baboon  experiments  we  con¬ 
sidered  peak  sled  acceleration  and  stopping  distance.  (Those  experi¬ 
ments  were  carried  out  on  the  Daisy  Decelerator.)  We  found  no  sta¬ 
tistical  evidence  that  the  occurrence  of  fatal  injury  depended  on 
stopping  distance.  Thus,  only  peak  sled  acceleration  was  used  to  pre¬ 
dict  injury.  In  the  baboon  experiments,  each  animal's  head  was  taped 
prior  to  experimentation,  precluding  the  use  of  initial  head  conditions 
as  predictor  variables. 

A  model  was  fit  to  the  combined  data  from  both  sets  of  experiments, 
and  the  baboons  shown  to  be  significantly  more  susceptible  to  the  effects 
of  -Gx  acceleration.  The  difference  in  LD^s  for  the  two  species  was 
estimated  to  be  about  24  g.  However,  the  differential  effect  of  accel¬ 
eration  was  found  to  be  the  same  for  these  species.  In  terms  of  toler¬ 
ance  limits,  this  implies  that  the  mean  tolerance  differs  for  the  two 
species  but  that  the  within-species  variability  in  tolerance  is  the 
same.  This  supports  an  hypothesis  advanced  earlier  by  Smith,  and  if 
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